What do we learn from the local geometry of glass-forming liquids? 
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We examine the local geometry of a simulated glass- forming polymer melt. Using the Voronoi 
construction, we find that the distributions of Voronoi volume P(vy) and asphericity P(a) appear 
to be universal properties of dense liquids, supporting the use of packing approaches to understand 
liquid properties. We also calculate the average free volume (vf) along a path of constant density 
and find that {v/) extrapolates to zero at the same temperature To that the extrapolated relaxation 
time diverges. We relate (v/) to the Debye- Waller factor. 



A fundamental mystery in the formation of glasses is 
the relationship of liquid structure to dynamics. While 
it has long been appreciated that many equilibrium and 
transport properties of dense fluids depend on the space 
available for molecular motion, the lack of methods to 
accurately compute or measure "free volume" and other 
measures of local structure has limited the development 
of this perspective for understanding fluid properties. 
Rahman Q suggested that the Voronoi cells (where a 
cell is defined as the sub-volume whose interior is closer 
to a specific atomic or molecular vertex than to any other 
vertex) may provide useful information about the local 
molecular environment. Furthermore, Voronoi volume 
vy has been considered a possible measure of the lo- 
cal free volume that might be correlated with liquid re- 
laxation (|, ||, and there has been some tentative evi- 
dence to support this proposition j|, [5) . Unfortunately, 
an unequivocal definition of free volume has been elu- 
sive, making it difficult to quantitatively test these ideas. 
However, for hard spheres, free volume can be rigor- 
ously defined and quantitatively related to the equation 
of state §0,|. 

In this Letter, we focus on the local geometry of a sim- 
ulated glass-forming polymer melt using a Voronoi anal- 
ysis and a free volume approach based on a mapping to 
hard spheres. From the Voronoi approach, we find the 
striking result that the distribution of cell volumes P(vy) 
and the distribution of cell asphericities P(a) are univer- 
sal functions of temperature T and density p over a wide 
range of the liquid state, independent of the interaction 
potential. While we are unable to quantitatively explain 
the origin of this universality, the existence of regularity 
in the structure of liquids is likely connected to the suc- 
cess of liquid state theories that focus on packing effects 
arising from core repulsion Utilizing a free volume 
definition related to hard spheres, we find that the ex- 
trapolated average free volume (w/) appears to vanish 
at the same temperature To where the extrapolated re- 
laxation time r of the coherent intermediate scattering 
function diverges. By relating (vf) to the Debye- Waller 



(DW) factor (u 2 ), we find that r ~ exp(u 2 /{u 2 )) (u is a 
constant), as experimentally observed [|10[ , and predicted 
by several models jll| . 

Our results are based primarily on molecular dynam- 
ics simulations of a melt containing 100 chains of "bead- 
spring" polymers, each chain consisting of 20 monomers. 
All monomers interact via a force-shifted Lennard Jones 
potential, truncated at 2.5, in standard reduced units. 
Neighboring monomers along a chain also interact via 
a FENE spring potential. The FENE parameters are 
k = 30 and i?o = 1-5, and are chosen to avoid crys- 
tallization. The dynamics of this model (with only the 
potential shifted) was studied in ref. p2| . Most of our 
simulations are in the range 0.35 < T < 4.0 at constant 
density p = 1.0, and as a result the average Voronoi cell 
volume (vy) = 1-0 is constant. We also consider several 
additional densities to demonstrate the generality of our 
results. For reference, we fit r to the power law form 
r ~ (T/T MCT — l) -7 expected from mode coupling the- 
ory (MCT) [[13| , and estimate the crossover temperature 
Tmct = 0.35±0.01, where the uncertainty represents the 
range of T for which a power law fits well to the data. 
A. Voronoi Cell Volume and Asphericity 

The sensitivity of the Voronoi analysis to local struc- 
ture is particularly helpful for understanding changes 
along an isochoric cooling path, since the changes in lo- 
cal environment are more subtle than along an isobaric 
path. We implement an efficient algorithm for generating 
the Voronoi tessalation and the Delauney simplices fu) . 
From these calculations we first consider the statistical 
properties of the Voronoi cell volumes vy. This leads to 
the striking result that the distribution of Voronoi cell 
volumes P(yv) appears to be universal, where the stan- 
dard deviation a 2 = (v v ) — (vy) 2 decreases upon cooling 
(Fig. |]). The scaling of P(vy) suggests the existence of 
a single underlying geometrical structure of the system, 
and that system specifics, such as temperature, density, 
and interaction potential, are absorbed into the average 
and variance of the distribution. To further test this pos- 
sibility, we calculate P(vy) for p = 0.9,0.95,1.05, and 
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FIG. 1: Distribution of Voronoi cell volumes, shifted by the 
average voronoi volume (iv) and scaled by the root- mean- 
square deviation a v . The figure includes data from 21 differ- 
ent state points of variable T and p (ranging from a very high 
T liquid to supercooled states), as well as data from 7 state 
points of simulated liquid water, and one state point of sim- 
ulated silicon. Most of the spread in the collapse is from the 
water data, for which we have inferior statistics. The inset 
shows the temperature dependence of a v along the p = 1.0 
iso chore. 

1.1 at T = 1.0 and find that P(vy) for different densities 
also scales to the same master curve. As an even more 
stringent test, we calculate P{vv) for configurations of 
simulated liquid water Jl5[ and silicon; both liquids are 
anomalous due to the presence of directional bonding in 
the liquid state. We find that P{vy) for water and sili- 
con collapses to the same universal curve as the polymer 
data, supporting the existence of a single underlying dis- 
tribution of Voronoi volumes for dense liquids. 

Such universal behavior suggests that a limit theo- 
rem |lfj| governs P(v\r) for a wide range of possible in- 
teractions. The distribution is best fit by a log-normal 
distribution, but the deviations of the fit exceed the qual- 
ity of the data collapse. This universal packing behavior 
likely explains why models based on short-range packing 
repulsion can successfully account for many liquid state 
properties ||. This may serve as a point of departure 
for theoretical investigation of the origin of the scaling, 
which is beyond the scope of the present paper. We also 
point out that scaling must break down at sufficiently low 
density or high temperature, since particles are randomly 
distributed in the ideal gas limit. 

The scaling parameter <j v follows a power law a v ~ 
T o.23±o.oi ( inget of p ig ^ Since av ig the fi uc t U ation 

in Voronoi volume, and compressibility kt is a measure 
of density fluctuations, it is tempting to relate the quan- 
tities Q. However, a v measures volume fluctuations 
the size scale of single particles, too small to expect a 
fluctuation relation of the form kt — ^/(wyfcflT) to 
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FIG. 2: Distribution of Voronoi cell asphericities, shifted 
by the average asphericity (a) and scaled by the root-mean- 
square deviation o a - The figure includes data from the same 
state points as Fig. [j]. The inset shows the T dependence of 
(a) along the p = 1.0 isochore. 

hold jL8j . On the other hand, a v should be related to 
the local restoring force of the fluid felt by a particle, 
since this controls the susceptibility of the fluid to local 
fluctuations. The restoring force is quantified by Qq, the 
first non-trivial coefficient in an expansion of the velocity 
auto-correlation function |js| |. Indeed, it has been found 
that f2rj scales approximately as T 1 / 4 pj| , which we also 
confirmed for our system. This suggests a connection be- 
tween <j v and vibrational dynamics (given the uncertain- 
ties in the exponent values). At fixed T = 1, we further 
observe a v ~ p~ 3,4 over the narrow range 0.9 < p < 1.1. 
In the cases of simulated water and silicon, our data is 
not able to confirm or exclude power-law scaling of a v . 

We further characterize local geometry by the aspheric- 
ity a of the cells, which we define by the ratio of the 
radius of a sphere with the same volume as the cell to 
the distance between the Voronoi vertex and the nearest 
cell face. If the cell is spherical, the ratio is one, while 
for any non-spherical cell, the ratio is greater than one. 
On cooling, a decreases by « 3% over the range of our 
simulations, and thus Voronoi cells become slightly more 
spherical on average (inset of Fig. [2]), as would be ex- 
pected if the packing becomes more regular on cooling. 
The distribution of asphericities P(a), like P(vy), also 
scales onto a single master curve using the same scaling 
method (Fig |2|) . The width of the distribution, quantified 
by the standard deviation a a , monotonically decreases 
with decreasing T (not shown), also expected if packing 
becomes more regular on cooling. We find that a a follows 
a power law whose exponent changes slightly at T ss 1, 
roughly the "onset" temperature for slow dynamics p(J 
and spatially heterogeneous monomer motion plj . How- 
ever, the lack of corroborating evidence from the simu- 
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lations of water and silicon suggests that the result for 
the polymer system is merely coincidental. Similar cal- 
culations for other simple and molecular liquids will be 
useful to identify the common features of cell geometry 
in liquids. 

B. Free Volume and Dynamics 

We next focus our attention on possible connections 
between local structure and dynamics by examining the 
free volume. There have been numerous definitions of 
free volume |^, 0], and thus it is not clear what the 
most appropriate and/or useful definition is. However, 
in the case of hard spheres, vt of a particle can be un- 
ambiguously defined as the volume over which the center 
of a sphere can translate, given that all other spheres in 
the system are fixed (see, e.g. Fig. 1 of ref. J22j, where 
an algorithm for the calculation can also be found) . This 
"rattle" free volume is rigorously related to the equation 
of state for hard spheres, as well as the void volume 
and thus is an attractive definition. 

In the case of soft-core repulsion, such a definition of 
Vf is not straightforward, since the average distance of 
closest approach depends on T. As done in ref. |2J| we 
define a T dependent effective hard sphere diameter, or 
Boltzmann diameter ub, by U(r = as) = E c , where U 
is the pair potential and E c is the average kinetic en- 
ergy of collision of an isolated pair of particles. Since E c 
for this system is not known, we estimate E c = 2ksT, 
the exact result for hard spheres. Additionally, we use 
only the non-bonded (Lennard-Jones) part of the poten- 
tial to determine <jb, since the bonding term has little 
effect on the distance of closest approach. By using the 
T dependent definition of <tb, it is possible that (vf) van- 
ishes along paths of constant density, such as we study 
here, in contrast with free volume definitions tied to the 
macroscopic density. 

We first consider the T dependence of the average free 
volume per monomer {vf) and r (Fig. ||). As expected, 
(vf) decreases on cooling while r increases. More sig- 
nificantly, we find that (i>f) may be well approximated 
by a power law (v f ) ~ (T/T - where T = 0.22 
and r\ = 1.46, if we allow both parameters to be free. For 
reasons that we will discuss, we expect n = 3/2, which re- 
sults in a best fit value of To = 0.20. We fit r to the Vogel- 
Fulcher-Tammann (VFT) equation r ~ expL4/(T - T )] 
and independently obtain an extrapolated divergence 
temperature To = 0.20, which is typically slightly below 
the laboratory defined T g . Hence, r seems to diverge at 
the same T at which (vf) extrapolates to zero, consistent 
with the possibility that the glass transition is related 
to vanishing of free volume, a central tenet of most free 
volume approaches pij. By eliminating the T depen- 
dence from the expressions for (vf) and r, we obtain the 
parametric relation 

T~eM(v/{v f }) 2/3 } (1) 
where v is a constant. This relation is similar to, but dis- 
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FIG. 3: (a) The relaxation time r of the intermediate scatter- 
ing function F(q,t); we define r at the time at which F(q,t) 
has decayed by a factor of e _1 at the q- value corresponding 
to the first maximum in S(q). The inset shows the quality of 
the fit to the VFT form, (b) Average free volume per particle 
(vf). The inset is a log-log plot to emphasize the power law 
behavior 



tinct from the Doolittle expression r ~ exp(v / v f) , later 
rationalized by the Cohen- Turnbull models that assume 
an exponential distribution of free volume Q . 

In order to explain the power-law dependence of (vf) 
on T, and hence the unexpected relation between r and 
(vf), as well as provide an experimental connection, we 
must better understand the physical origin of (vt). Since 
the free volume measures the space over which a parti- 
cle can move before encountering the exclusion volume of 
neighboring particles (in the hard sphere mapping), it is 
natural to expect that free volume and the Debye- Waller 
factor, a measure of cage size, might be related. To test 
this, we define the DW factor (u 2 ) = (r 2 {t = 1.022)), 
where t — 1.022 is the approximate time of the crossover 
from ballistic to caged motion of the mean-squared dis- 
placement (r 2 (t)}. In Fig. || we make a parametric plot of 
(v-f) and (w 2 ) 3 / 2 and find a linear proportionality, sup- 
porting the hypothesis that free volume should be re- 
lated to the DW factor. As a more stringent test, we 
consider the distributions P(vf) and P(u 3 ), where P(u 3 ) 
is calculated in the same fashion as the van Hove cor- 
relation function. The inset of Fig. |^ shows that both 
distributions are nearly exponential, but that stretching 
occurs at large volume. Exponential decay of P(vf) at 
large Vf is essential to recover the Doolittle relation in 
the Cohen- Turnbull formulation of free volume theory; 
since the Doolittle relation does not appear to hold for 
our data, the significance of the deviation from exponen- 
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FIG. 4: Parametric plot of the Debye- Waller factor (u 2 ) raised 
the the 3/2 power (such that is has units of volume) as a 
function (Vf). The line indicates a least-square fit. The inset 
shows the distributions P(u 3 ) and P(vf). The data for P{u 3 ) 
is arbitrarily shifted to emphasize the similarity of the distri- 
butions. The bold line in an exponential exp(n/0.2 3 ), plotted 
as a guide to the eye. 



tial decay is unclear. More importantly, the similarity of 
P(vf) and P(it 3 ) for sufficiently large volume allows us 
to interpret v/ asa measure of the DW factor, and apply 
ideas that have already been developed relating (u 2 ) and 
bulk relaxation properties. 

We substitute (w 2 ) 3 / 2 for (vf) in Eq. |]and obtain 



exp(u 2 /(u 2 )). 



(2) 



Several arguments have been put forth for this phe- 
nomenologically observed relation jl0| , based on the idea 
that the effective force constant localizing a particle is in- 
versely proportional to (u 2 ) and directly proportional to 
the energy barrier height O]. The T dependence of (u 2 ) 
is controlled by two factors: (i) inertial energy, result- 
ing in a term proportional to ksT, counterbalanced by 
(ii) a nearly T independent elastic energy proportional to 
Good 3 , where Goo is the shear modulus and <5 3 ~ 0{vy)- 
This immediately yields (vf) ~ (T/T — 1) 3 ^ 2 , where 
T oc GoaS 3 /k B , since (v/) oc (u 2 ) 3 / 2 . Such a relationship 
is consistent with the VFT expression, as demonstrated 
by simple substitution. Additionally, the Adam-Gibbs 
expression r ~ exp(A/T5 con f), successfully applied to a 
variety of model liquids [ p5[ , suggests a non-trivial rela- 
tionship to the configurational entropy S'conf ; specifically, 
P^conf OC (u 2 ) or TSconf oc (w/) 2/3 . This offers an area of 
future consideration. 

We have focused our attention on a limited subset of 
the predictions of free volume theories. We did not find 
support for several additional hypotheses. Specifically, 
we also find that (i) free volume percolates at T well 
above T g and (ii) that that no significant correlation ex- 
ists between local volume of a specific monomer and its 



mobility. More details on these results will be provided 
in a future publication. 

We thank J. Baschnagel, W. Kob, F. Sciortino, and V. 
Novikov for helpful discussions. 
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